Magnetic field induced finite size eff'ect in type-II superconductors 
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We explore the occurrence of a magnetic field induced finite size effect on the specific heat and 



^ ; 

' correlation lengths of anisotropic type-11 superconductors near the zero field transition tempera- 

, ture Tc. Since near the zero field transition thermal fluctuations are expected to dominate and 

^ ^ ' with increasing fleld strength these fluctuations become one dimensional, whereupon the effect of 

^ I fluctuations increases, it appears unavoidable to account for thermal fluctuations. Invoking the scal- 

. ing theory of critical phenomena it is shown that the specific heat data of nearly optimally doped 
YBa2Cu307_5 are inconsistent with the traditional mean-field and lowest Landau level predictions 

, of a continuous superconductor to normal state transition along an upper critical field Hc2 (T). On 

' ■ the contrary, we observe agreement with a magnetic field induced finite size effect, whereupon even 

' the correlation length longitudinal to the applied field H cannot grow beyond the limiting magnetic 

, length Lh oc ^J(^o/H. It arises because with increasing magnetic field the density of vortex lines 

O ' becomes greater, but this cannot continue indefinitely. Lh is then roughly set on the proximity of 

O I vortex lines by the overlapping of their cores. Thus, the shift and the rounding of the specific heat 

5^ , peak in an applied field is traced back to a magnetic field induced finite size effect in the correlation 

^L(' length longitudinal to the applied field. 

^3 ; 

. The superconductor to normal state transition in conventional low Tc materials appears to be well described by the 
' Ginzburg-Landau mean-field approximation. Because of the large correlation volume in these materials, the region in 
which critical fluctuations are important is too small to be accessible experimentally. In contrast, with the discovery 
I of superconductivity in the cuprates, a new era started [1]. Indeed, marked deviations from mean-field behavior have 
. been observed over a temperature range of the order of 10 K above and below Tc [2-22]. Theoretical expectations of 
' the kind of critical behavior which might be observed are: (i) If fluctuations in the vector potential can be ignored, 
then the zero-fleld transition belongs to the universality class of the three dimensional XY-model, as is the superfluid 
transition in "^He. In an applied magnetic field the critical behavior is then equivalent to that of uniformly rotating 
''He near the superfiuid transition [22,23]; (ii) When fluctuations in the vector potential are included the charge 
of the Cooper pairs, entering via the effective dimensionless charge e — ^/A = 1/k, charged critical behavior is 
. expected to occur in which both the correlation length ^ and the magnetic penetration depth A grow with the same 
' critical exponent by approaching Tc from below [24-33]. However, in extreme type-II superconductors where k >> 1 
. the effective charge e is very small. As a consequence the region close to Tc, where the system crosses over to the 

■ regime of charged fluctuations, becomes too narrow to access. For instance, optimally doped YBa2Cu307_5, while 
, possessing an extended regime of critical fluctuations, is too strongly type-II to observe charged critical fluctuations 

■ [2-22]. In strongly type-II superconductors (k >> 1) the crossover upon approaching Tc is thus initially to the 
"Y^! • critical regime of a weakly charged superfluid where the fluctuations of the order parameter are essentially those of an 

d ' uncharged superfluid or XY-model. Furthermore, there is the inhomogeneity induced flnite size effect which renders 
^ the asymptotic critical regime unattainable [20,21]. However, underdoped cuprates appear to open a window onto the 
I ' charged critical regime because k becomes rather small in this doping regime. Here the cuprates undergo a quantum 
superconductor to insulator transition in the underdoped limit [16,22] and correspond to a 2D disordered bosonic 
system with long-range coulomb interactions. Close to this quantum transition Tc, Xat and ^ab scale as Tc (x 
O ■ ^ iab [16,22], yielding with the dynamic critical exponent z — I [16,22,34-36], Kab oc a/T^- Recent measurements of 
the magnetic in-plane penetration depth of underdoped YBa2Cu3 06.59 clearly uncovered critical behavior associated 
with a charged critical point, in which both the coherence length and the magnetic penetration depth grow by 
approaching Tc from below with the same critical exponent [37] . Thus, as far as static zero field critical phenomena 
are concerned, there is little doubt that near optimum doping the observable critical behavior of bulk YBa2Cu307_5 
is governed by the three dimensional (3D) XY universality class. Accordingly, we expect that the critical behavior in 
an applied magnetic field is equivalent to that of uniformly rotating ^He near the superfluid transition [22,23]. The 
singular part of the free energy per unit volume should then scale as [16,22] 



abj 



for a magnetic field He applied parallel to the c-axis. G (z) is a universal scaling function, ± = sgn (t) = 
sgn {T /Tc — 1), the zero-field in-plane correlation length with critical amplitude C^jg, 7 = 'C^fco/^rf) anisotropy, 
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and is a universal constant, fixed by the 3D-XY universality class. The fluctuation contribution to the magneti- 
zation m = —dfs/dHc scales then as 

while the singular part of the specific heat, c{Hc, T) = c {He, T) p/ks — —d'^fs/d\t\^ adopts the scaling form 

c{He,T) = — \t\-''F^{z), ^^(^"fco) =-Q±a{l-a){2-a) = {R^)'\ 3iy = 2 - a, (3) 
a 7 

p denotes the density and c{Hc,T) is in units of {erg/ {gK)). In the 3D-XY universality class are the universal 
quantities Q^, R^, v and a given by [38] 

~ 0.36, R- ~ 0.82, ~ 1.06, ~ 2.28, z/ = (2 - a) /3 ~ 0.67 (5) 

In the presence of a sufiiciently small magnetic field He, the specific heat is then expected to have a singular part 
which exhibits the scaling behavior 



where 



(6) 



The magnetization data [5,10] and the zero-field specific heat measurements of YBa2Cu307_5 [4,9,16] agree well 
with these predictions. In a nonzero applied field, one can test the scaling form (6) of the specific heat by the extent 
to which data for c{Hc,T) Hc^"^^ collapse to a common curve when plotted as a function of x. Here, matters are 
complicated by the fact that a different kind of scaling behavior 

Z{He, T) = R {xl) , XL = (7) 

is expected when only the lowest Landau level (L) is significantly occupied [39,40]. Here T — Tc2 {He), or equivalently, 
He = Hc2 {T) is the upper critical field of the Ginzburg-Landau theory. Since a in Eq.(6) is very small, and 2f « 4/3, 
the two predictions are rather hard to distinguish [41]. Some authors argue that lowest-Landau- level scaling works 
just as well as critical-point scaling [42-52]. 

Theoretically, the scaling form (6) is an unambiguous prediction of the theory of critical phenomena and ought to 
be observed sufficiently close to the zero- field critical point. On the other hand, lowest-Landau-Ievel scaling, relics on 
the assumption that the correlation length longitudinal to the applied magnetic field, diverges along the line Tc2 {H) 
[53], whereupon a continuous phase transition from the superconducting to the normal state is predicted to occur. 
Accordingly, the behavior of the correlation length longitudinal to the applied field is essential to verify the lowest 
Landau level prediction. In this context it is instructive to rewrite the scaling variable z (Eq.(l)) in the form 

z = ^ = Lh, = V^o/ {aH,), (8) 

with a ~ 3.12 [20], related to the average distance between vortex lines [20,22,54 56]. The scaling function G {z) is 
then identical to that of a system with finite extent Lh^ in the a6-plane [57,58]. As a consequence, fluctuations which 
are transverse to the applied field are stifle and the corresponding correlation lengths cannot grow beyond 



- V*o/ {aHk) =LH„i^j^k. (9) 

Hence, the fiuctuations of a bulk superconductor in a magnetic field are longitudinal to the field and for this reason 
one dimensional, as noted by Lee and Shenoy [59]. Noting that fluctuations become more important with reduced 
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dimensionality, one expects that the remaining fluctuations, which are longitudinal to the applied field, remove the 
mean- field transition at Tc2 {He)- Indeed, thermal fluctuations destroy the ordered phase in one dimensional systems 
with short-range interactions. Furthermore, calculations treating these interactions within the Hartree approximation 
[39,40,54,60], and generalizations thereof [41,61,62], find that the correlation length longitudinal to the applied field 
remains bounded as well. In this case the correlation length ^~ {Hc,t) adopts the scaling form 

^,{H,,t)=^f,\trS^{x), x= (^2^£2Lj \t\. (10) 

The scaling function must behave as 

5=^ (a; = ±oo) = 1, 5=^ (x ^ ±0) = sq l^r , (11) 



so that for He 0, {Hc,t) oc \t\ " and for t ^ 0, {Hc,t) oc y^^o/Hc ■ Thus, the divergence of the correlation 
length ^~ {He, t) is removed and it adopts a maximum at Tp {He) < Tc, yielding the line 



H,c{t)= ^ ;_ , \tr. (12) 

On the other hand, the specific heat scales according to Eq.(3) as 

Z{H,,t) = ^\t\-" F^{x). (13) 
In this case the scaling function behaves as 

F± {x = ±00) = 1, F± {x = ±0) = /± , (14) 

so that for He — ^ 0, c{Hc,t) cx and for t ^ 0, c{Hc,t) cx ($o/-f^c)~"^^''- Noting that the ratio 

{x) I {F^ {x))"^" tends to constant values in the hmits a; ^ and x ±00 the relation between the {t, He) 
and the absolute value of the specific heat \c{Hc,T)\, 



{Hc,t) = 



^ ,~ , .~.,vja 'S'^ (2;) 



a 



c{HeM''^ (15) 

^ ^' (F±(x))'^/" ^ ' 



as obtained from Eqs.(lO) and (13), reduces in these limits to 

€e(i?c,i)oc|?(ife,i)r/". (16) 

Thus, when this scenario holds true, the specific heat probes essentially the correlation length longitudinal to the 
applied field. As a remnant of the singularity at in zero field it should exhibit a so called finite size effect [57,58], 
resulting in a smooth peak around Tp because the correlation length (t, H^ cannot grow beyond (tp, (Eq.(12)). 
Furthermore, given experimental data for c{Hc,t), this scenario can be verified. Indeed in the limits He and 
t ^ ±0 the critical behavior \c{He,t)\'"^°' oc 1*1"'' with Co/Co - 2-28 (Eqs.(5)) should hold. This allows to 
circumvent the aforementioned difficulties associated with the comparison of scaling functions with the prediction 
of the lowest-Landau-level approach. Indeed, if there is a magnetic field induced finite size effect on the correlation 
length longitudinal to the applied field, the transition is rounded and the assumption of an upper critical field He2 is 
not justified. 

Here wc analyze the specific heat data of Roulin et al. [49] to verify the magnetic field induced finite size scenario. 
In Fig. la we depicted the data for the YBa2Cu307_5 single crystal in terms of c{He,t) cx c{He,T) /T — B vs. 
t = T/Te-l with B = 0.1717 J/K^gat and = 92.77 K. The solid and dashed lines area {He, t) (x c {He, T) /T-B = 
with A- = -0.0672 J/K2gat, A+ = -0.0685 J/(K2gat), and a = -0.01 (Eq.(5)). Using the relations 
A- (J/(K gat)) = A-Tea, A- (cm^^) = W"^ / {ksYgat) A' (J/(K gat)), and Vgat = 8 cm^ we obtain for A' the 
estimate A^ — 5.64 10^*^ cm~^ and with the universal relations (3) and (5) for the correlation volume and the critical 
amplitudes of the zero field correlation lengths the estimates 

(Go) ' ^co ^ 978 A3, :^ 18.43 A, ~ 2.88 A, (17) 
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using 7 = 6.4. As a remnant of the zero-field singularity, there is for fixed field strength a smeared peak adopting its 
maximum at Tp which is located below Tc- As Tp approaches Tc, the peak becomes sharper with decreasing He and 
evolves smoothly to the zero-field behavior, smeared by the inhomogeneity induced finite size effect, arising from the 
limited extent Lab,c of the homogeneous domains along the a6-plane and c-axes. Since Tp decreases systematically 
with reduced field down to He — 0.25 T, corresponding to Lh^ < \/^o/ iflHc) = 512 A, the magnetic field sets at and 
above He = 0.25 T the limiting length. On the contrary, when Lh^ > Lab, the inhomogeneities set the limiting length. 

In this case, tp — {£,abo/-^abY^'^ , independent of the applied field. Thus, the field dependence of tp is a characteristic 
feature of a magnetic field induced limiting length. The line Hep (t) is shown in Fig. lb. The solid line is Eq.(12) with 

*o/ (Cbo)') = 210 T, yielding with Cw - 18-43 A (Eq.(17)), 

Xp~2.21, (18) 
and Zp = Xp^^^" ~ 0.55, which agrees well with the previous estimate Zp ~ a^^^'^ = 3.12^^^^ ~ 0.57 [20]. 




FIG. 1. a) c{Hc,t) (X c{Hc,T) /T - B vs. t for YBa2Cu307-i derived from Roulin et al. [49] with B = 0.1717 J/K^gat 
and Tc = 92.77 K; The solid and dashed lines are c{Hc,t) oc c(Hc,T) /T - B = \t\~" with A" = -0.0671 J/K^gat, 
A+ = -0.0684 J/K^gat, and a = -0.01. b) Hep (t) vs. t . The solid line is Hep (t) = 210 Itf" with u = 0.67, corresponding to 
Eq.(12) with $o/ {xl- iCmf) = 210 (T). 

In Fig. 2 we displayed the scaling plot c{He,t) oc {c{He,T) /T — B) vs. t/H^^'^^ derived from the data 
of Roulin et al. [49] shown in Fig.la. Noting that according to Eq.(13), {c{He,T) /T - B) cx c{He,t) = 
A^F^ (x) /a, where x cx t/H^^'^, this plot uncovers essentially the scaling function (x), whereupon the data 
points should fall on a single curve, as they apparently do, when plotted versus x cx t/Hc^'^'^ . The solid and dashed 
curves indicate the asymptotic behavior in the limits x cx t/ Hc^'^" — s- ±0 (Eq.(14)), while the arrow marks tp/ h}/'^^ — 
, 2 \ ^1"^" 

Xp I (^afco) /'^o ) ' where the peak in the specific heat adopts its maximum value. However, as aforementioned, it 
is difficult to distinguish different models on the basis of such scaling functions. 
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FIG. 2. c{Hc,t)\t\~°' oc {c{Hc,T) /T ~ B)\t\" vs. t/HV^" derived from the data of Roulin et al. [49] shown in 
Fig. la. The solid and dashed lines are —0.0654 (^—t/Hl^'^"^ and —0.0658 (^—t/Hl''^'^^ , respectively. The arrow marks 
tp/H^^^'' ~ -0.0184. 

However, in view of Eqs.(15) and (16), giving the relationship between the fluctuation contribution to the specific 
heat and the correlation length longitudinal to the applied field, {t, He), the magnetic field induced finite size scenario 
can be verified by considering the plot \c{Hc, t)\'^^°' cx |c (T, He) /T — B\"^°' vs. t, predicted to probe {t, He). In Fig. 3 
we show |c(i?cji)r^" cx: \c{T,Hc) /T — B['^°' vs. t derived from the data of Roulin et al. [49] shown in Fig. la. The 
solid and dashed line mark the leading zero field critical behavior in terms of |c(i?c,i)r^" oc \c{T,Hc) /T ~ B\'^°' — 
\ty^ , with r~/r+ = 2.28, consistent with the universal ratio ^^/CJo — 2.28 of the 3D-XY universality class 
(Eq.(5)). This confirms that in the scaling regime considered here \c{Hc,t)\'^'^ probes essentially the correlation 
length longitudinal to the applied magnetic field, {t,Hc), so that Eq.(16) applies. The rounded peak in zero field 
reveals then an inhomogeneity induced finite size effect, while the smeared peak in finite fields, its shift and broadening 
with increasing field strength discloses the magnetic field induced finite size effect on the correlation length longitudinal 
to the applied field. Because {t,Hc) o: \c(Hc,t)\'^^" cx \c{T,Hc) /T — B\^^°' the field dependence of tp, where the 
correlation length adopts its maximum value set by the magnetic field, coincides with the tp {He), where the specific 
heat reaches its maximum value. In this context it is important to recognize that the magnetic field dependence of tp 
is a unique consequence of the magnetic field induced finite size effect. Indeed, when inhomogeneities set the limiting 

length then tp = {^abal ■^"■'^^^^ ^ independent of the applied field. 

When the magnetic field induced finite size effect scenario is correct, the occurrence of the effect is not restricted to 
temperatures below T^. It is particularly dramatic at T^, where in a homogeneous system in zero field the correlation 
lengths are infinite. In an applied field the scaling form (10) yields the prediction 



where we used the definition (8 for the limiting magnetic length Lh^- In Fig. 3b we displayed ^c{Hc,t = 0) oc 
\c{Hc,t — 0)1''^" cx \c{Hc,Tc) /Tc — B\'^^°' vs. He- The consistency with the solid line, which is Eq.(19) in the form 
ic{H^,t = Q) cx |c(ifc,i = 0)|''/"°'|c(iJc,rc)/rc-B|''/" = 3.210^9i^c~^/^ reveals again that an applied magnetic 
field leads to a finite size effect in the correlation length longitudinal to the field. 




(19) 
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FIG. 3. a) |?(i?c,t)r''" oc \c{H,,T) /T ~ Bl"'" % 
The solid and dashed lines are \c{Hc,T) /T ~ B\'"'" 



b) \5{H,,t = 0)\""^ (X \c{H,,T,)/T,-B 
\c{H,,n)/n~B\-^°' =3.2 10^9//, 



-1/2 



«. t derived from the data of Roulin et al. [49] shown in Fig. la. 
= r^ltp" with P- = 4.7 10''** , r+ = r^/2.28 , and u = 0.67. 
He derived from the data shown in Fig. 3a. The solid line is 



For fields applied parallel to the a-axis, the transverse correlation lengths and are according to Eq.(8) bounded 
by ^j£,a£,^ — Lh^ ■ When the magnetic field induced finite size scenario holds true, the correlation length longitudinal 
to the apphed field, {Ha,t), should be bounded as well. In analogy to Eq.(lO) the longitudinal correlation length 
adopts then the scaling form 



-l/2i/ 



\t\. 



(20) 



with the limiting behavior given in Eq.(ll). Thus, the divergence of the correlation length is removed and (Ha,t) 
adopts at Tp (Ha) < Tc a, maximum, yielding the line 



Hpa {t) 

Furthermore, in analogy to Eq.(16) the relation 



f<,(ii„ocx|H(i/„i)r/", 



(21) 



(22) 



between the longitudinal correlation length and the specific heat should hold in the limits a; ^ and x ±oo. 
In Fig.4a we depicted \c{Hab,t)\'"^" cx \c{Hab,T) /T- B]"/" vs. t derived from the data of Roulin et al. [49]. The 
solid and dashed line mark the leading zero field critical behavior in terms of \c{Hc^ t)\^^°' cx |c (T, He) /T — Sj''^" = 
r='= 1^1"'', with r^/r+ — 2.28, consistent with the universal ratio ^abol^lbo — ^-^^ 3D-XY universality class 

(Eq.(5)). This consistency confirms again that in the scaling regime considered here |c (ffob, t)]"^" probes the corre- 
lation length ^af, (t, Hab) longitudinal to the applied field. The rounded peak in zero field reveals an inhomogeneity 
induced finite size effect, while the smeared peak in finite fields, its shift and broadening with increasing field strength 
disclose the magnetic field induced finite size effect in ^ab {t, He). Indeed, from Fig. 4b it is seen that the field depen- 
dence of tp where the correlation length S^ab (t, He) adopts its maximum value, set by the magnetic field, is consistent 

with tp = -O.OOSiHll^" which is Eq.(21) with ((a;2'^CfcoCo) /*o)^^^" = 0.0034 (T), resulting from the C^j, given by 
Eq.(17) and Xp = 2.21 (Eq.(18). 
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FIG. 4. a) IciHab,^^" oc \c{Hab,T) /T - B]"^" vs. t derived from the data of Roulin et al. [49] with B = 0.1717 J/K^gat 
and Tc = 92.77 K. The solid and dashed lines are \c{Hah,T) /T - B\"^" = \t\'" with T" = 5.05 10™ , r+ = r^/2.28 , and 

u = 0.67. b) vs. Hab- The sohd line is tp = -O.OOMHH'^" , which corresponds to Eq.(21) with {{xp" ^^^^Q^) /$o)^^^'^ = 0.0034 
T. 



To summarize our result for an anisotropic type-II superconductor, we have shown that near the zero field transition 
temperature superconductivity is in a magnetic field subjected to a field induced finite size effect. The crucial ingredient 
for a finite size effect is an energy gap in the excitation spectrum of fluctuations. In the present case it is the discrete 
set of Landau levels. Indeed, there is the formal analogy with the Landau levels of a charged particle moving in circular 
orbits in the plane perpendicular to the applied field at the cyclotron frequency. As a consequence, the fluctuations 
which are transverse to the field are stiff and have a length scale Lh oc y/^o/H. Hence, the fluctuations of a bulk 
type-II superconductor become one dimensional and are longitudinal to the applied field, as noted by Lee and Shenoy 
[59] . Because fluctuations become more important with reduced dimensionality, one expects then that the interaction 
of these fluctuations remove the mean-field transition at Tc2 (He), because thermal fluctuations destroy the ordered 
phase in one dimensional systems with short-range interactions. The absence of this transition is further supported 
by calculations treating the fluctuations within the Hartree approximation [39,40,54,60], and generalizations thereof 
[41,61,62]. They suggest that the correlation length longitudinal to the appHed held remains bounded as well. Invoking 
the scaling theory of critical phenomena we confirmed this prediction. We have shown that the specific heat data 
of Roulin et al. [51] clearly reveals a magnetic field induced finite size effect in the correlation length longitudinal to 
the applied field. Accordingly, there is no evidence for a phase transition line Tc2 (H) near the zero field transition 
temperature T^. 
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